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ABSTRACT 

Silicate features arising from material around pre-main sequence stars are useful probes of the star 
and planet formation process. In order to investigate possible connections between dust processing and 
disk properties, 8-13 /im spectra of 34 young stars, exhibiting a range of circumstellar environments 
and including spectral types A to M, were obtained using the Long Wavelength Spectrometer at 
the W. M. Keck Observatory. The broad 9.7 fim amorphous silicate (Si-0 stretching) feature which 
dominates this wavelength regime evolves from absorption in young, embedded sources, to emission 
in optically revealed stars, and to complete absence in older "debris" disk systems for both low- 
and intermediate-mass stars. This is similar to the evolutionary pattern seen in ISO observations of 
high/intermediate-mass YSOs. The peak wavelength and FWHM are centered about 9.7 /mi and ^2.3 
/im, respectively, corresponding to amorphous olivine, with a larger spread in FWHM for embedded 
sources and in peak wavelength for disks. In a few of our objects that have been previously identified 
as class I low-mass YSOs, the observed silicate feature is more complex, with absorption near 9.5 /tm 
and emission peaking around 10 /tm. Although most of the emission spectra show broad classical 
features attributed to amorphous silicates, small variations in the shape/strength may be linked to 
dust processing, including grain growth and/or silicate crystallization. For some of the Herbig Ae 
stars in the sample, the broad emission feature has an additional bump near 11.3 /tm, similar to 
the emission from crystalline forsterite seen in comets and the debris disk /3 Pictoris. Only one of 
the low-mass stars, Hen 3-600A, and one Herbig Ae star, HD 179218 clearly show strong, narrow 
emission near 11.3 /im. We study quantitatively the evidence for evolutionary trends in the 8-13 /im 
spectra through a variety of spectral shape diagnostics. Based on the lack of correlation between 
these diagnostics and broad-band infrared luminosity characteristics for silicate emission sources we 
conclude that although spectral signatures of dust processing are present, they can not be connected 
clearly to disk evolutionary stage (for optically thick disks) or optical depth (for optically thin disks). 
The diagnostics of silicate absorption features (other than the central wavelength of the feature), 
however, are tightly correlated with optical depth and thus do not probe silicate grain properties. 

Subject headings: circumstellar matter — stars: formation — stars: pre-main-sequence — planetary 
systems: protoplanetary disks — infrared: stars 
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1. INTRODUCTION: SILICATES AND STAR FORMATION 

Circumstellar disks are found around most, if not all, 
young stars. At least 5 % of these circumst ellar disks 
are likely planet-forming l)Marcv et al.l 12000) and it is 
commonly believed that our solar system formed from a 
solar nebula similar to circumstellar disks found around 
young stars (Russell et al., Met. and Early Solar Sys- 
tem II in press). Initially these disks are composed of 
gas and dust from the interstellar medium. Theoret- 
ical models predict that collisions between grains and 
subsequent coagulation result in modified grain size dis- 
tributions and eventually the growth of planetesimals 
ijWeidenschilling fc Cuzzi|ll993j) . Howev er, observational 
evide nce of this process remains elusive ijBeckwith et alJ 
2000j . One way to gain insight into the process of grain 
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growth, planctesimal evolution, and eventually planet 
formation, is to constrain the effects of these processes 
on the observed mid-infrared spectra from circumstellar 
dust. Grain growth and changes in composition due to 
dust processing in circumstellar disks can have dramatic 
effects on the observed spectrum of the dust in emis- 
sion and absorption. In this contribution, we attempt to 
quantify those effects and place them in the context of 
modern theories of disk evolution and planet formation. 

Dust in the interstellar medium, and thus in the early 
stages of star formation, consists primarily of oxygen- 
rich silicates and c arbon-rich grains (for summary, see 
iPollack et al.l lL994|) from a wide range of sources, in- 
cluding stellar atmospheres, novae, supernovae. These 
silicates are either olivines (Mg 2:c Fe( 1 _ a .)Si04), ranging 
from fayalite (x — 0) to forsterite (x = 1), or pyrox- 
enes (Mg^Fe^-^SiOs), ranging from ferrosillite (x = 0) 
to enstatite (x = 1). Silicate condensates formed in su- 
pernovae and during star formation are primarily Fe-rich 
and amorphous, but at high temperatures, Mg-rich crys- 
talline silicates (i.e., pure forsterite and enstatite) are 
produced. The fraction of pyroxenes versus olivines de- 
pends on both temperature and pressure. For a given 
pressure, the pyroxenes (e.g., enstatite) dominate at low 
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temperatures and becomes less abundant with increas- 
ing temperature (at equilibrium). The olivines (e.g., 
forsterite) form a significant fraction of the mixture only 
for a narrow region in pressure-temperature space. 

The low temperatures of the ISM are expected to result 
in amorphous mixtures of olivines and pyroxenes with 
olivine the dominant component. Indeed, strong smooth 
features near 9.7 (Si-0 stretching mode) and 18.5 /im 
(O-Si-0 bending mode) are observed in absorption along 
lines of sight through the diffuse interste llar medium 
ijWhittet et alJll997t iBouwman et alJl200I|) . Molecular 
clouds also possess a mixture of olivines and pyroxenes, 
which are expected to remain amorphous due to cosmic 
ray impacts; a broad structureless band centered near 
10 /im is observed toward several o bjects (j Wfilner et^rt.1 
1981 IPeeourie fc Papoularl [l985t IWhittet et alJ 119881 
Bowev et al.lll998|K Observations bv lBowev fc Adamsonl 
l)2002[) suggest a shift from amorphous olivine dominated 
grains in diffuse-medium environments to amorphous py- 
roxenes in molecular clouds and YSOs, indicating that 
some dust processing is occurring during the early stages 
of star formation. Further processing is evidenced by 
the presence of crystalline silicates in the dust surround- 
ing some main sequence stars, with amorphous olivine 
and forsterite as t he dominant components (see e.g., 
iKnacke et al] 1*1993). Comets also contain bo th amor- 
phou s and crystalline silicates (see review by iWoodenl 
2002), with crystalline fractions that are too high (up to 
30% by mass) to be accounted for by the ISM (w here 
the crystalline fraction is < 2%: lKemner et al.l*2004ft . in- 
dicating that crystallization occurred in the proto-solar 
nebula. 

Differences in the silicate spectra of interstellar ver- 
sus solar-system dust suggest an evolutionary sequence 
in which the grain composition changes from amorphous 
olivine to crystalline forsterite during the star and so- 
lar system formation process. Dust grains in accretion 
disks around young stars therefore must be significantly 
modified with respect to those of the ISM and are thus 
excellent laboratories for the study of in situ crystalliza- 
tion of silicates. If crystallization proceeds via annealing, 
then silicates in the outer disk are expected to accrete 
icy mantles and retain cores that are similar to ISM dust 
(i.e., amorphous olivine and pyroxene), while in the in- 
ner disk (T > 800 K), microcrystalline lattice structures 
on 0.001 /i m size scales can form within silicate grains 
( Gail 1998). If annealing continues for a sufficiently long 
time, or occurs in a sufficiently hot environment, crystal- 
lization spreads from these critical nuclei over the entire 
grain. At disk radii corresponding to temperatures above 
~ 100 K, 0.1-1 ttm grains should be completely crystal- 
lized l)Gai]l"L998') and larger dust grains should possess 
regi ons of localized crys tallization detectable in IR spec- 
tra ll.Taeger et a.l .1 1*1994). Increases in temperature are 
also predicted to result in conversion from pyroxenes to 
olivines. Thus we would expect both crystallization and 
composition (i.e., the content of Fe vs. Mg and olivine vs. 
pyroxene) of the silicates in disks to be dependent on ex- 
posure to stellar radiation (e.g., location in the disk, disk 
geometry, stellar type) and/or the duration of processing 
(e.g., amount of disk turbulence, disk age). 

High spectral resolution studies of the warm dust emis- 
sion near 10 /im are vital in understanding the com- 
position and processing of dust in the planet-forming 



region (1-10 AU) of circumstellar disks around young 
stars. The dominant feature in the 8-13 /im spectral 
region accessible from the ground is from amorphous 
olivine (peaking at 9.5 /im), with weaker features aris- 
ing from silica (SiC>2; 8.6 /im) and crystalline silicates 
including enstatite (MgSiC>3; 9.2, 10.4 /im) and forsterite 
(Mg2Si04; 10.2, 11.3 /im) 6 . The transition from amor- 
phous olivine to crystalline forsterite results in a pre- 
dicted shift of the silicate emission feature from 9.7 to 
11.3 /im. However, analysis of the 10 /im silicate fea- 
ture is complicated by the fact that the strength and 
shape of the emission can not be uniquely correlated 
with abundance or composition. The presence of organ- 
ics on grain surfaces results in increased opacity at A < 
6 /tni and decreased contrast between t he silicate band 
at 10 /im and the adjacent continuum ifD'Alcssi o et alJ 
l200l . Gr ain size has a similar effect and for large 
grains (a max > 10 (jm) the silicate band s disappear al- 
most entirely (Wolf & Hillenbrand 2003). Additionally, 
iFabian et alJ l)2001l) suggest that, grain shape also needs 
to be taken into account in the interpretation of sili- 
cate spectra, particularly for strong crystalline bands, for 
which a combination of spherical and ellipsoidal grains 
may be required. Disk morphology and optical depth 
also influence the shape /strength of the silicate feature 
l)Meeus et al.ll200lt IBouwman et al.l *2003). Despite diffi- 
culties in constructing unique interpretations, a number 
of ground- and space-based studies of the disks around 
T Tauri and Herbig Ae/Be (hereafter referred to as 
HAEBE) stars have been conducted, yielding important 
constraints on the evolution of circumstellar material. 

Early studies of young low-mass stars llCoherJ Il980t 
[ Cohen fc WittehornllT985l IWoodenl 1*19941 IHanner et all 
I1995L Il998[) show a single, rather sharp peak be- 
tween 9 and 10 /im, resembling molecular clouds and 
the diffuse ISM . ISO spectra of young low-mass stars 
in Chameleon l|Natta et alJ 120001) display 10 /zm sili- 
cate emission that is significantly broader and stronger 
than that from the diffuse ISM, but consistent with 
the emission predicted by a simple radiation trans- 
fer model of <1 /im amorphous pyroxene and/or 
olivin e grains in a flared disk. More recent surveys 
re.gJSitko et a.l.l1999U200(TUSvlyester fc Manning J200(t 
Hond a et all I20 03: Mecus et al. 2003: P rzvgodda et alJ 
I20031 Ivan Boekel et al.ll2003|) . however, have found ev- 
idence in a few low-mass sources for a more trapezoidal 
shaped 8-13 /im amorphous silicate feature or even a pro- 
nounced 11.3 /tm crystalline silicate (forsterite) bump, 
indicating that at least some dust processing has taken 
place in these objects. 

HAEBE stars were studied more extensively with ISO 
and, based on 2-45 /im ISO spectra of HAEBE stars 
spanning a rang e of ages, an evolutionary sequence has 
been suggested ijMeeus fc Waelken£*ll999j) . In this se- 
quence, rising featureless dense cloud core spectra pro- 
ceed to strong amorphous 10 /im silicate absorption from 
embedded protostars. These are followed by optically 
thick protoplanetary disks with spectra showing warm 
10 /im silicate dust emission and prominent cold dust 
emission longward of 15 /im, then spectra with only 
warm dust interpreted as arising from disks undergoing 

6 The peak wav elengths described above are for 0.1 fim grains 
IMeeus et al.ll200ll) . 
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planet formation, and finally the most evolved sources 
with nearly photosphcric spectra, indicating disk dis- 
persal. In addition to this general evolutionary trend, 
ISO detected crystalline silicate emission bands (11.3, 
27.5, 33.5, 35.8 and 70 /Ltm) tow ard a modest sam- 
ple of HAEB E star+disk systems ijMolster et alJ 119991: 
M eeus et all 1200 1|) . The longer wavelength emission 
bands are highly dependent on the structure of the dust. 
These studies show large variations in dust crystallinity 
and compositi on among YSOs of similar age and stel- 
lar p roperties ijvan den Ancker et al.ll2000t IMeeus et al.l 
120011) . indicating the necessity of large samples to test 
scenarios of disk evolution. 

Objects that may be in transition between pri- 
mordial dust and debris dust (generated through 
collisions of planetesimals), or which are fully de- 
bris, have also been investig ated l|Svlvester et alJll996t 
iSvlvester fe Mannings! I200"o() revealing generally photo- 
spheric or continuum excess spectra. The 10 um sp ec- 
trum of the 12±\ Myr old l|Zuckerman et all l200l|) A 
star (3 Pictoris, however, exhibits a remarkable resem- 
blance to that of solar system grains and comets, sug- 
gest ing^jJ L j]iixture_ofj3rystalline and amorp hous ma- 
terial ijKnacke etHTI Il993t iReach et all l2003|) . Direct 
imaging of (3 Pictoris shows a warped disk that is po- 
tentially caused by collisions of planete simals and may 
be sh aped by recent planet formation llKalas fc Jewittl 
119951 iMouillet et aflll997l: lAueereau et al.ll200H) . Spa- 
tially resolved 10 um spectroscopy of [3 Pictoris 
l)Weinberger et aLl l2"003') revealed that the silicate emis- 
sion arises only from dust within 20 AU of the star, 
compared to the 77 AU emitting region in the contin- 
uum. Furthermore, the mixture of amorphous and crys- 
talline dust appears to remain constant, suggesting that 
the crystallization process may be acting uniformly as a 
function of radius within the planet formation region. 

Observations of silicate emission toward young stellar 
objects in a variety of evolutionary stages can provide 
valuable information about the modification of dust dur- 
ing solar system formation. In the present study, 8-13 
/im spectra were obtained for 34 low- and intermediate- 
mass young stars, ranging from embedded protostars to 
debris disk systems. Additional photometry at 10.7 /im 
was obtained for a subset of the sample in order to facili- 
tate flux calibration of the spectra. In §2, the sample se- 
lection criteria are discussed. The spectroscopic and pho- 
tometric observations and data reduction are described 
in §3. SEDs compiled from the literature are presented 
in §4, along with implications regarding the evolutionary 
categories within the sample. In §5, the 8-13 /im spectra 
are discussed qualitatively and in §6 a quantitative anal- 
ysis of the silicate emission spectra is presented. Con- 
clusions and suggestions of future work are discussed in 
§7. The appendices contain the references for the SEDs 
(A) and a short discussion of variability in the silicate 
feature (B). 

2. SAMPLE SELECTION 

Spectra in the 10 /jm region of a moderate sample of 
optically visible HAEBE, T Tauri, as well as more heavily 
embedded YSOs were acquired. The goal was to use the 
spectra to characterize the composition of the dust, par- 
ticularly the degree of crystallinity of the silicate grains 
present. Objects with varying age and mass were chosen 



in order to put this information in context with the stellar 
evolutionary process for both low- and high-mass stars. 
The sample includes 7 lo w-mass class I objects chosen 
from IKenvon et all l)1993j) and 3 well known high- mass 
YSOs: the Becklin Neugebauer object (BN), NGC 2024 
IRS2, and Mon R2 IRS3. These observations probe the 
composition of largely unprocessed grains and absorp- 
tion from amorphous silicates is expected. There are ten 
optically visible class II T Tauri and HAEBE stars in the 
sample, most of which have very well characterized SEDs. 
ISO observations indicate that at least two HAEBE stars 
in this intermediate-age sample, HD 163296 and HD 
179218 , possess detectable amoun t s of crystall i ne sil - 
icates l|van den Ancker et all 120001: IMeeus et al"ll2001|) . 
Seven older main sequence stars with debris disks were 
also included, which possess optically thin IR continua 
and ages older than 6 Myr. These are the most likely 
candidates for detection of crystalline silicates, provided 
warm dust is present. Our source list and observing dates 
are provided in Tabled 

3. PHOTOMETRIC AND SPECTROSCOPIC 
OBSERVATIONS 

8-13 /im spectra/photometry were obtained using the 
Long Wavelength Spectrometer (LWS) at the W. M. 
Keck Observatory between 1999 August and 2000 De- 
cember (see Table LWS provides diffraction- limited 
(FWHM = 0."22) imaging (10" field of view, 0."08 
pixel" 1 ) and long slit (7") spectroscopic (R = 100-1400) 
capabilities in the 3-25 /im wavelength range. Photom- 
etry was obtained with the 10.7 /im filter (AA = 10.0- 
11.4 /zm) for sources observed during the nights of 2000 
Feb 20-21 and December 9, under poor and variable 
seeing conditions: 0."3-0."6 at 10 /im. The Nwide fil- 
ter and low-resolution grating (LRES) were used to ob- 
tain 8-13 /im spectra, with R = 100, at a dispersion 
of 0.037 /jm pixel -1 . Slit widths of 3 and 6 pixels were 
used, resulting in 0."24 and 0."48 apertures, respectively. 
Sources were imaged prior to each spectroscopic obser- 
vation to ensure optimal placement within the slit. For 
most spectroscopic observations, additional calibration 
scans were obtained using the Keck routine LSEC, in 
which the data is chopped between an ambient black- 
body source (T « 274 K) and the sky. Darks were taken 
at the end of each night with exposure times equal to 
those on source. The raw "six-dimensional" LWS data, 
consisting of images for two chop-nod pairs, were coad- 
ded into two-dimensional images, using the IDL routine 
LWSCOADD 7 . 

The IRAF/PHOT task was used for the photomet- 
ric data reduction. Apertures of radius 0."96 (12 pix) 
were used for the photometry and the sky background 
was measured in 2. "0-3. "2 radius (25-40 pix) annuli. 
Flat-fielding was found to increase the scatter in the pho- 
tometry and therefore not performed. A curve-of-growth 
correction of -0.12 ± 0.01 mag (based on standard star 
measurements) to a 1".92 (24 pix diameter) aperture was 
applied to each star, resulting in magnitudes within ^5% 
of the infinite aperture value. Standard star observations 
were used to obtain extinction curves for each night. We 

7 LWSCOADD was written by Gregory D. 
Wirth and is provided by the observatory at 
http://www2.kcck.hawaii.edu/inst/lws/lwscoadd.html 
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TABLE 1 
Observations 









RA 


DEC 


Observation 


Reference 


Aperture 




Source 


(J2000) 


(J2000) 


Date 


Star 


(pixels) 


LSEC? 


49 Cet 




01 


34 


37.78 


-1 5 40 34 9 


2000 Nov 8 


gl Tau 


g 




HD 17925 




02 


52 


32.129 


1 9 4fi 10 Q79 


9000 Foh 9i 


HR 1 01 7 

111 v. 1U1 I 




y 


HD 22049 




03 


32 


55.844 


OQ 97 9Q 7zL4 


9000 npr- i o 


UR 1 01 7 

llll. 1U1 ( 


u 


y 


IRAS 040164 


-2610 


04 


04 


42.85 


4-26 18 56.3 


2000 Dec 9 


HR 617 


6 


y 


IRAS 041084 


-2803B a 


04 


13 


52.9 


4-28 11 23 


2000 Dec 10 


HR 1017 


6 


V 


IRAS 041694 


-2702 


04 


19 


59.24 


4-27 09 58.6 


2000 Dec 10 


HR 1017 


6 


V 
J 


IRAS 041814 


-2654A/B a 


04 


21 


11.42 


4-27 01 08.9 


2000 Dec 10 


HR 531 


6 


V 
J 


IRAS 042394 


-2436 


04 


26 


57.1 


4-24 43 36 


2000 Dec 10 


HR 1017 


6 


V 


IRAS 042484 


-2612 


04 


27 


56.7 


4-26 19 20 


2000 Dec 9 


HR 617 


6 


v 

J 


IRAS 042644 


-2433 a 


04 


29 


07.68 


4-24 43 50.1 


2000 Dec 9 


HR 1017 


6 


v 


Haro 6-10A 




04 


29 


24.39 


4-24 33 02.1 


2000 Dec 10 


HR 1017 


6 


V 


Haro 6-10B 




04 


29 


24.39 


4-24 33 02.1 


2000 Dec 10 


HR 1017 


6 


V 


IRAS 042874 


-1801 


04 


31 


33.6 


4-18 08 15 


2000 Dec 9 


HR 1708 


6 


V 


IRAS 042954 


-2251 


04 


32 


32.07 


4-22 57 30.3 


2000 Dec 9 


HR 617 


6 


v 


AA Tau 




04 


34 


55.2 


4-24 28 52 


2000 Feb 21 


HR 2990 


6 


V 
J 


IRAS 043254-2402 b 


04 


35 


33.0 


4-9ZL 08 1 A 


9000 Dec Q 








LkCa 15 a 




04 


39 


17.8 


4-99 91 03 


1 QQQ Nnv 30 9000 Nrw £ 


A r 


6 3 


n, 


IRAS 043814 


-2540 a 


04 


41 


12.48 


4-25 46 37 1 


2000 Dec 10 


HR 531 


g 


y 


IRAS 044894-3042 a 


04 


52 


06.9 


+30 47 17 


9nnn PpK 91 


TTR 1 70K 
nix j. i uo 


g 


y 


GM Aur a 




04 


55 


10.2 


-u^n 91 


9000 Nmr R 

AUUU 1\UV o 


q Tau 


3 




MWC 480 




04 


58 


46.27 


4-2Q 'SO ^7 


iqqq M nv ■in 


a Tau 


u 




BN 




05 


35 


14.17 


-05 22 23 1 


2000 Dec 10 




fi 


y 


NGC 2024 IRS2 


05 


41 


45.8 


-01 54 30 


2000 Dec 10 


HR 2943 


6 


y 


Mon R2 IRS3 


06 


07 


47.8 


-06 22 55 


2000 Dec 10 


HR 2943 


6 


y 


HD 233517 




08 


22 


46.71 


4-53 04 49.2 


2000 Nov 8 


a Aur 


3 


ii 


Hen 3-600A 




11 


10 


27.9 


-37 31 52 


2000 Feb 21 


HR 4786 


6 


y 


HD 102647 




11 


19 


03.578 


4-14 34 19.417 


2000 Feb 19 


HR 4534 


6 


y 


HR 4796A 




12 


36 


01.032 


-39 52 10.219 


2000 Feb 21 


HR 4786 


6 


y 


IRAS 14050-4109 a 


14 


05 


05.7 


-41 09 40 


2000 Feb 21 


HR 4786 


6 


y 


HD 163296 




17 


56 


21.29 


-21 57 21.9 


1999 Aug 23, 2000 Jun 20 


P Oph 


6,3 


n,y 


HD 179218 




19 


11 


11.254 


4-15 47 15.630 


2000 Feb 21 


HR 5908 


6 


y 


WW Vul a 




19 


25 


58.75 


4-21 12 31.3 


2000 Jun 20 


e Cyg 


3 


y 


HD 184761 




19 


34 


58.97 


4-27 13 31.2 


2000 Nov 8 


a Tau 


3 


ii 


HD 216803 




22 


56 


24.0529 


-31 33 56.042 


2000 Dec 9 


HR 1017 


6 


y 



a Residuals of incomplete subtraction of the 9.7 /im telluric ozone feature are present in the spectrum. 
b Photometric observations only. A 8-13 /an spectrum was not acquired for this source. 



derived a typical night-to-night variation in the atmo- 
spheric zero-points of ^20% and an average atmospheric 
extinction correction of -0.32 ± 0.07 mag airmass" 1 . The 
root-mean-square scatter in the photometry of the stan- 
dards is 0.06 mag for 2000 February 20 and 21 and 0.13 
mag for 2000 December 9. 

The spectroscopic data reduction was performed using 
the NOAO IRAF/TWODSPEC and ONEDSPEC 
packages. After applying pixel masks and flat fields, the 
images were then divided by those produced from the 
LSEC calibration scan obtained closest in time to the 
observation, to correct for variations in the spectral re- 
sponse of the grating. The detector's dark current was 
found to be negligible over the integration times used 
and therefore not subtracted from the data. Spectral 
extraction was performed via standard methods using 
the IRAF task apall. For the purposes of wavelength 
calibration and removal of the telluric ozone absorption 
feature at ~9.5 /im, each extracted spectrum was aligned 
with and divided by a standard star spectrum (see Ta- 
bic Q] for details) using the IRAF telluric task. In some 
cases (as noted in TableQ]), variations in seeing and atmo- 
spheric transmission versus wavelength resulted in small 
differences in the shape of the telluric ozone absorption 
feature observed in the target and calibrator spectra, and 
residuals can be seen in the resulting divisions. After 



telluric corrections were performed, the resulting spec- 
tra were multiplied by a blackbody emission spectrum 
appropriate to each calibrator to remove any induced 
slope in the 8-13 /im region. Absolute flux calibration 
was performed using photometric observations at simi- 
lar wavelengths culled from the literature or obtained in 
this study (as indicated in Tabic [2J , with scale factors 
obtained from the spectrum integrated over the photo- 
metric bandwidth. Flux calibrated spectra are shown in 
Figures ITMH1 The spectra are overplotted onto the SEDs 
to indicate the success of this flux calibration method. 

4. SPECTRAL ENERGY DISTRIBUTIONS (SEDS) 

Continuum fluxes for the sources in our sample were 
collected from the literature in order to form the spec- 
tral energy distributions shown in Figures HTBI The data 
were not de- reddened, but, as the Ay for all but the 
class I sources is small (< 3 mag), this should have a 
modest effect on fluxes for wavelengths longer than 1.2 
/im (J-band). In Figures 13161 the SEDs for stars with 
known spectral types are overlaid with emission arising 
from a blackbody at the stellar temperature found in 
the literature for each source (see Table |3J) by scaling 
the peak of the blackbody curve to match the SED at 
that wavelength. More quantitatively, stellar luminosi- 
ties were found by calculating a color excess and A„ for 
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TABLE 2 
Flux calibration of spectra 





A 


Flux 


AA 






Source 


((im) 


(Jy) 


(A™ 1 ) 


Instrument 


Reference 


49 Cet 


10.8 


0.20 ± 0.04 


8.0-13.6 


KeckII,OSCIR 


1 


HD 17925 


10.7 


0.92 ± 0.09 


10.0-11.4 


Keck,LWS 


2 


HD 22049 


12 


9.5 ± 0.4 


9.15-14.85 


IRAS 


3 


IRAS 04016+2610 a 


10.7 


2.8 ± 0.2 


10.0-11.4 


Keck,LWS 


1 




10.3 


2.5 ± 0.1 


9.65-10.95 


IRTF 


5 


IRAS 04108+2803B a 


10.7 


1 + 3 


10.0-11.4 


Keck,LWS 


1 




10.3 


0.57 ± 0.09 


9.65-10.95 


IRTF 


5 


IRAS 04169+2702 


12.0 


0.75 ± 0.05 


9.15-14.85 


IRAS 


3 


IRAS 04181+2654A/B 


12.0 


0.36 ± 0.05 


9.15-14.85 


IRAS 


3 


IRAS 04239+2436 


12.0 


1.71 ± 0.09 


9.15-14.85 


IRAS 


3 


IRAS 04248+2612 


10.7 


0.230 ± 0.009 


10.0-11.4 


Keck,LWS 


1 


IRAS 04264+2433 


10.7 


0.41 ± 0.02 


10.0-11.4 


Keck,LWS 


1 


HARO 6-10A 


10.0 


4.3 b 






6 


HARO 6-10B 


10.3 


10 ± 2 


9.65-10.95 


IRTF 


5 


IRAS 04287+1801 


10.7 


5.0 ± 0.4 


10.0-11.4 


Keck,LWS 


4 


IRAS 04295+2251 


10.7 


0.69 ± 0.03 


10.0-11.4 


Keck,LWS 


4 


AA Tau 


10.7 


0.44 ± 0.03 


10.0-11.4 


Keck,LWS 


2 


IRAS 04325+2402 c 


10.7 


0.060 ± 0.003 


10.0-11.4 


Keck,LWS 


4 


LkCa 15 


9.6 


0.50 ± 0.15 


9.5-9.7 


ISO,SWS 


7 


IRAS 04381+2540 


10.1 


0.23 ± 0.05 


9.65-10.95 


IRTF 


5 


IRAS 04489+3042 


10.7 


0.31 ± 0.02 


10.0-11.4 


Keck,LWS 


1 


GM Aur 


10.1 


0.5 


7.55-12.65 


IRTF 


8 


MWC 480 


9.6 


8.7 ± 2.6 


9.5-9.7 


ISO,SWS 


7 


BN 


10.7 


88 


10.0-11.4 


Keck,LWS 


9 


NGC 2024 IRS2 


10.5 


19 ±1 


8.0-13.0 


MANIAC, ESO 


10 


Mon R2 IRS3 


10.0 


165 ±5 


9.65-10.95 


IRTF 


11 


HD 233517 


10.1 


0.439 


9.595-10.505 


UKIRT.Berkcam 


12 


Hen 3-600A 


10.7 


0.73 ± 0.06 


10.0-11.4 


Keck,LWS 


2 


HD 102647 


10.7 


5.51 ± 0.41 


10.0-11.4 


Keck,LWS 


2 


HR 4796A 


10.7 


0.22 ± 0.02 


10.0-11.4 


Keck,LWS 


2 


IRAS 14050-4109 


10.7 


0.23 ± 0.02 


10.0-11.4 


Keck,LWS 


2 


HD 163296 


9.69 


18.6 ± 0.5 


8.84-10.54 


lm,ESO 


13 


HD 179218 


10.7 


17 ± 2 


10.15-11.25 


O'brien, bolometer 


14 


WW Vul 


9.6 


2.3 ± 0.5 


9.5-9.7 


ISO,SWS 


7 


HD 184761 


12 


2.41 


9.15-14.85 


IRAS 


3 


HD 216803 


10.7 


1.53 ± 0.082 


10.0-11.4 


Keck,LWS 


2 



Re ferences. — (1) IJayaward hana et aT| l200H ; (2) IMetchev et alJ IgO Oa (3*1 IBe ichman et alJ 
19881 (4) this paper; (5) IMyers et al.| 1 19871 (6) extrapolation from ILeinert fc Maaj|T989J; (7) 
Thi et alJ l200l (8) IKenvon et alj|1993t (9) Hillenbrand; unp ublishe d, CCD; (10) IWalsh et alJ 
200 U: rilt IKoresko et alJ I1993t (12) ISkinner et all 119951 (13) IBerrilli et all 119921; (14) 
Lawrence et al.H1990l 

a Photometric errors are larger in the data acquired in this study, therefore 10.3 fim photometry 
from (5) were used for flux calibration. 
b Flux for Haro 6-10B is an estimate by Lcincrt & Haas (1989). 

c Photometric observations only. A 8—13 /im spectrum was not acquired for this source. 



the given spectral type, applying a bolometric correction 
to the V or I magnitude (as indicated), and converting 
to bolometric luminosity using the given distance (see 
Table |3 for details). For comparison, stellar luminosities 
from the literature are also presented. These are similar 
to the luminosities derived here; a notable exception is 
HD 179218, for which the V-band magnitudes reported 
in the literature vary by up to 1.3 magnitudes. For con- 
sistency, we will use the luminosities derived here, but 
note that there may be errors for stars with large vari- 
ability. 

The IR luminosity and the ratio L/^/L* (indicative of 
dust geometry for optically thick disks and of midplane 
optical depth for optically thin disks) were estimated 
from the SEDs collected for our sample in the follow- 
ing manner. First, the wavelength (X on set) at which the 
photometric flux is greater than twice that of the stellar 



blackbody was evaluated 8 . The infrared luminosity of 
the disk Lm was then calculated by integrating the SED 
between X on set and 1 cm and subtracting the integral of 
the stellar blackbody over the same region. The resulting 
IR-to-stellar luminosity ratios (L/fl/L*), as well as the 
wavelengths of the onset of excess radiation (X on set), are 
presented in Table As the spectral types and Ay for 
the embedded YSOs are uncertain, only total luminosi- 
ties were calculated; resulting values of L tot are presented 
in Table H 

The SEDs were used to establish the continuum around 
the 8-13 /im silicate emission feature, which is extremely 
important for analysis of the feature shape and the 

8 A onse t is calculated to give the reader an impression of where 
the SED begins to deviate from photospheric and was chosen to 
be the point at which the photometric flux is greater than twice 
that of the stellar blackbody as this minimizes the errors due to 
variability in the photometric measurements. 
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TABLE 3 

Stellar parameters from the literature and calculations from SEDs: revealed disks 





Distance 


Spectral 


L* 


T eff 


Age 








L*.boi d 


\ c 

s^onset 




spectral 


Source 


(pc) 


Type 


(L0) 


(K) 


(Myr) 


Ref. a 


Method b 


A c 


(L©) 


(H 




index g 


49 Cet 


U± 


A1V 


23 44 




7 R 


1,2 


V/B-V 


21 




zo 




1 94 


HD 17925 


10.4 


K1V 




5000 


80 


3 


V/B-V 





0.43 


60 


<6.0(-5) 


2.64 


HD 22049 


3.2 


K2V 




5000 


330 


3 


V/B-V 





0.33 


60 


9.7(-5) 


2.61 


AA Tau 


140 


K7 


0.71 


4060 


2.1 


2 4 


I/V-I 


1.0 


0.60 


2.2 


0.48 


0.84 


LkCa 15 


140 


K5:V 


0.74 


4350 


2.0 


2,5 


I/V-I 


1.1 


0.86 


1.2 


0.49 


0.90 


GM Aur 


140 


K7 


0.724 


4060 


1.8 


2,4 


J/J-H h 


0.06 


0.73 


6.9 


0.43 


0.48 


MWC 480 1 


131-140 


A3e 


11-32 


8460-8890 


4.6 


6,7 


V/B-V 


0.25 


14 


1.6 


0.29 


0.53 


HD 233517 


23 


K2 


0.20 


4500 




8 


I/V-I 


0.30 


0.12 


10 


0.12 


0.62 


Hen 3-600A 


50.0 


M4 




3290 


10.0 


9 


I/V-I 


0.10 


0.43 


11 


9.1(-3) 


0.62 


HD 102647 


11.1 


A3V 




8590 


240 


3 


V/B-V 





14 


60 


<3.9(-6) 


2.59 


HR 4796A 


76 


AOV 


35 


10000 


8.0 


10,11 


I/V-I 





26 


25 


4.6(-3) 


0.92 


IRAS 14050-4109 


140 


K5 




4405 




12 


I/V-I 


0.81 


0.78 


12 


8.6(-2) J 


0.68 


HD 163296 


122+H 


AlVe 


30.2 


9332 


4.0 


6 


V/B-V 


0.34 


32 


1.2 


0.28 


0.58 


HD 179218 


240±I« 


B9e 


316 


10471 


0.1 


6 


V/B-V 


0.59 


110 


2.2 


0.32 


0.13 


WW Vul 


550 


A3e/B9V 


43 


8600 




14 


V/B-V 


0.95 


30 


1.6 


0.50 


0.57 


HD 184761 


65 


A8V 




7500 




15 


V/B-V 


0.09 


7.4 


2.2 


0.10 


1.38 


HD 216803 


7.70 


K5Ve 




4555 


200 


16,17 


V/B-V 


0.06 


0.20 


60 


<7.0(-4) 


2.67 



References. — (11 ICoulson et alj|1998t (2llThT ct al. 2001; (3) Habine ct alJI200U (4) E artmann c t alj|1998t (5) IStrom et al.H198St (61 
van den Ancker et alJll99SI: (7) ISimon et at] 120011 (8) ISMimer et aij|1995t (9) ITorres et alj 120031 (1 0) Uura et al.H1993t (11) IStautter et aU 
19951 (12) IGregorio-Hetem fe Heterd 120021 (13) [Maimings fc Sargentl 120001 (14) INatta et alj 120011 (15) IMiroshnichenko et al.H1999t (161 
Barrado v Navascues ct al. 1997, (17) Santos ct al. 2004. 

Note. — Parameters on the left are from the literature as noted; those on the right arc calculated here. 
a Where two references are listed, the second is for the age. 

b Method is in the form of A/B-C, where B-C is the color used to calculate the stellar A^ and A is the band to which the bolometric correction 
was applied. 

C A„ were calculated from the spectral type given and color excess E(B-C), as defined in the Method column. 

d Stellar bolometric luminosities were calculated by applying a bolometric correction to the magnitude of band A (defined in the Method 
column) and converting to luminosity using the given distance. 
e A nset is the wavelength at which the stellar blackbody and a polynomial fit to the SED diverge. 

f Lj r is calculated by integrating the SED between Xonset and 1 cm and subtracting the integral of the stellar blackbody over the same region. 
L* — L*,6oi- 

Log(\ b F\ ) — Log(\ a Fx ) 

s The spectral indices are calculated following the method of lKenvon fc Hartmannl 11995H . spectral index = Log(X b )-Log(X ) — ~' ^ a ~ ^'^ 

am for most sources, excepting IRAS 14050-4109 (0.8 fj,m) and HD 216803 (3.6 /^m); A;, = 25 am for all sources. 
h J- and H-band data was used to calculate the luminosity for GM Aur, as the I-band magnitude seemed anomalously high, 
'd = 131 pc and T e ff = 8710 were used for the luminosity calculations for MWC 480. 

J The SED for IRAS 14050-4109 could not be extrapolated beyond the last observed point (100/im) and therefore this luminosity ratio is a 
lower limit. 



silicate comp osition. To establish t he continuum for 
HAEBE stars. Ivan Boekel et alJ (|2003l) fit each SED with 
a blackbody for the star and several blackbodies of differ- 
ent temperatures for the dust emissio n and interpolated 
between the two. INatta et al.l l)2000t ) fit a powerlaw to 
the edges of the silicate emission feature to remove the 
continuum in a sample of T Tauri stars. Upon exam- 
ination of the overlays of the spectra on the SEDs in 
Figure EE1 we nn( i that for our sample the continuum is 
adequately represented by connecting the two endpoints 
of the spectrum (averaging over 0.20 /Ltm near 8.3 and 
12.2 fjtm). 

In general, we would expect the shape of the 8-13 //m 
feature to be related to the SED and to the relative disk 
versus stellar radiation. Specifically, if changes in the sil- 
icate feature are due to processing of the dust, through 
growth and/or crystallization, then the shape of the sili- 
cate emission feature should be related to the evolution- 
ary stage of the source. Because the spectral index is 
an indicator of the evolutionary class, we expect that 
it may be related to the shape of the 8-13 /jm spec- 
tra. In the following sections we explore these trends, 



first through qualitative discussion and then through 
quantitative characterization of the spectral strength and 
shape. 

5. DISCUSSION OF THE 10 fiM SPECTRA 

The flux calibrated 0.1-10 5 fj,m SEDs and 8-13 fxm 
spectra are shown in Figures 1 1 If il with the continuum 
subtracted silicate features shown in the right column 
in Figures The sources are placed into categories 

based on the shape of the SEDs and 8-13 /im spec- 
tra. In general, the debris disks are found to have lit- 
tle or no evidence of silicate emission, the T Tauri and 
HAEBE sources show strong to moderate silicate emis- 
sion and the embedded sources show silicate absorption, 
or a combination of absorption and emission. Possible 
evolutionary trends are shown in Figures and |S] for the 
low-mass and intermediate-mass stars, respectively. The 
trend shown in Figure |H1 is similar to that suggested for 
HAEBE stars bv lMeeus fc Waelkensl <fl999h . The evolu- 
tion of low-mass star spectra and SEDs does not appear 
to be significantly different from intermediate-mass stars, 
with embedded sources showing amorphous olivine in ab- 
sorption, followed by stages of dust processing (including 
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Fig. 1. — SEDs and 10 /im spectra for low-mass stars where silicates are in absorption. The left panel displays the SEDs for each source 
(see appendix A for references) in units of Log AF^, with the 10 /jm spectra included for reference. The insets show enlargements of the 
8-13 fim regions with linear scaling on both axes. Open symbols are data for which a single source within a binary or multiple can not be 
distinguished. The continuum is in most (but not all) cases adequately represented by connecting the two endpoints of the spectrum and 
normalized spectra are shown in the right panel in units of (F„-F c )/F c (see text for details). The dashed line depicts the continuum level. 



grain growth and/or silicate crystallization) resulting in 
variations in the shape of the emission feature, and finally 
photospheric spectra for debris disks. One possible dif- 
ference between the low- and intermediate-mass stars in 
our sample is that for the intermediate-mass stars, crys- 
talline silicates have been found in sources which possess 
large near- and far-IR excesses, but the one low-mass 
star toward which crystalline silicates have been found, 



Hen 3-600A (see also iHonda et ail 12003 lUchida et alJ 
l2004j) . has an SED with no near-IR excess and dimin- 
ishing far-IR flux. There have been far fewer detections 
of crystalline silicates toward low-mass, with respect to 
high-mass, stars reported in the literature as well. It is 
not at all clear, however, that the surveys to date are 
drawn from representative samples. 
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-continued 



5.1. Silicate Absorption Spectra 

The silicate feature is clearly seen in absorption in 9 
low-mass embedded sources and in all 3 high mass proto- 
stars in our sample (see Figures ^ and EJ). The spectra of 
IRAS 04287+1801 and IRAS 04381+2540 (and possibly 
Mon R2 IRS3) possess a shallow absorption feature near 
11.2 /zm in addition to the 9.7 fim absorption seen toward 
the other embedded sources. The additional absorptio n 
feature has rarely been seen (e.g., iBooeert et all 12004]) . 
and its origin may be related to grain mantle ice con- 



stituents and not silicates (Boogert, private communica- 
tion). In 5 other low-mass embedded objects the silicate 
feature appears with absorption superposed on emission 
between 9 and 12 //m. In the most extreme case, IRAS 
04264+2433, the spectrum shows almost pure silicate 
emission. These spectra may be explained by a source ge- 
ometry in which cold dust located along the line of sight 
is absorbing against silicate emission from warmer dust 
located near the protostar. Radiative transfer modeling 
is required to explore this scenario and constrain whether 
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-continued 



the cold dust is in fact associated with these sources or 
associated with remnant unbound material in the dense 
core. 

The absorption spectra exhibit a variety of morpholo- 
gies and have been ordered in Figure^roughly according 
the strength of the absorption feature. This may corre- 
spond to an evolutionary sequence from strong ISM-like 
absorption through superposed absorption+emission and 
will be discussed more quantitatively in §6.1. 

5.2. Silicate Emission Spectra 

The silicate feature is seen in emission toward 5 T Tauri 
stars and 5 HAEBE stars in our sample (see Figures |3] 
and QJ. 8 out of 10 of these stars possess significant 
IR luminosities (Lj/j/L* > 0.1), with the exception of 
IRAS 14050-4109 and Hen 3-600A (L /fl /L, ~ 10~ 2 , sim- 
ilar to debris disks). The shape of the emission feature 
varies from smooth single peaks near 9.7 /tm, consistent 
with amorphous olivine, to more complex spectra with 
additional peaks near 10.2 and 11.3 /mi, indicative of 



enstatite and forsterite components. This suggests that 
silicate processing is occurring within disks around opti- 
cally revealed T Tauri and HAEBE stars. 

8-13 /tm spectra similar to that of amorphous olivine 
were observed toward 4/5 T Tauri stars (GM Aur, AA 
Tau, IRAS 14050-4109 and LkCa 15). The emission 
features are wider in comparison to absorption features 
in the ISM. This could either be due to the addition 
of ^10.2-11.3 /xm emission from crystalline silicates or 
broadening of the amorphous emission feature due to in- 
creased grain sizes. The absence of a strong crystalline 
component in the emission spectra indicates that the 
wider feature is more likely related to grain growth (dis- 
cussed below) than to crystallization. Crystalline sili- 
cates appear to be absent even in the two sources with 
evidence of disk clearing, IRAS 14050-4109 (with Ljij/L* 
~ 10~ 2 ) and GM Aur (possessing a gap at small radii, 
iRice et al.ll2?)0l . 

Possible evidence for crystallization appears in the 
spectra observed toward 4/5 HAEBE stars (MWC 480, 
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MONR2 IRS3 




Fig. 2. — SEDs and 10 fim spectra for high mass stars where silicates are in absorption, otherwise the same as Figure HI 



TABLE 4 

Calculations from SEDs: embedded objects 





Distance 


Spectral 


Source 


Ltot (pc) 


Reference Index a 



1.03 
1.09 
1.29 
0.31 
0.44 
1.10 
0.49 
1.00 
0.71 
1.56 
1.78 
0.23 
1.31 
0.31 
1.25 
1.74 
2.77 



References. — (lV Hartrnann 2002; (2) Bcichman et alll988t 
(3) IM vers et al J [19871: (4) IGezari et al.H1998t ( 5 ) IHaisch et Id) 
120011 f6) IPreibisch et ai1l20021 

Note. — The luminosities for several close sources are those 
of the combined sources due to the large IRAS beam. 

a The spectral indices are calculated as in Table 151 \ a = 2.2 fim 
for all sources; Aj, = 25 (im for all low-mass stars and 20, 25, and 
39 /an for BN, NGC 2024 IRS2, and Mon R2 IRS3, respectively. 



HD 163296, HD 184761, and WW Vul). HD 163296, HD 
184761, and WW Vul each possess an emission feature 
characteristic of amorphous olivine as well as a small 
peak at 11.3 /*m, consistent with crystalline forsterite 
(though this peak is m uch weaker than that seen for 
Picto ris or HD 100546. iKnacke et aljfl993t IMeeus et al] 
12001 . MWC 480 possesses a peak near 10.5 /im, which 
is consistent with crystalline enstatite. The three disks 
with 11.2 fim peaks are young intermediate-mass stars 
(~2 M©) with SEDs depicting substantial IR emission 
from warm dust, indicating that the component produc- 
ing the feature must be present at a young evolution- 
ary stage. Analysis of the crystalline silicate content in 
these disks, however, may be complicated by the 11.2 /im 
PAH emission band. Clear identification of the 11.3 fim 
feature as arising from crystalline silicates can only be 
obtained if there is evidence of their longer wavelength 
features, or if the presence of PAHs can be ruled out. 
HD 163296 and WW Vul have been observed with ISO 
and no evide nce of PAH emission in the 6-12 /zm region 
was d etected ijMeeus et al.l200lllAcke fc van den Anckerl 
l2004|k longer wavelength crystalline silicate features were 
detected toward HD 163296. In the same study, PAH 
emission at 6.2 fj,m was detected toward MWC 480, but 
similar features at 7.6 and 8.5 were not. No previous 
observations of HD 184761 are available. 

Finally, one low-mass T Tauri star (Hen 3-600A) and 
one HAEBE star (HD 179218), 10 Myr and 6 Myr of 
age, respectively, each show a strong peak at 11.3 /im 
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Fig. 3. — SEDs and 10 fim spectra for disks around low-mass T Tauri stars where the silicate 10 /an band is in emission. In the left 
panel, the SEDs are plotted in units of Log AF^ overlaid by dashed lines representing the stellar photosphere, as fit by the blackbody model 
described in the text. The inset shows an enlargement of the 8-13 fira region with linear scaling on both axes. The continuum is obtained 
by connecting the two endpoints of the spectrum. Continuum subtracted spectra in units of (F„-F c )/F c are shown in the right panel, with 
the continuum level depicted by dashed lines. 



that dominates the spectrum. The shape of the emission 
feature is consistent with that of crystalline forsterite (in- 
cluding the correlating small peak at 10.2 jim). These 
spectra are similar to that of the first widely recognized 
cryst alline silicate source, HD 100546 l|Bouwman et alJ 
2003). In contrast to the detection of crystalline sili- 
cates toward several HAEBE stars, to date crystalline 
silicates have been detected toward only a small number 



of -solar-mass s tars (TW Hya. iSitko et all [2000; VW 
Cha and Glass I, iMeeus et al.ll2003f) . Our sample seems 
consistent with this trend, with 1/5 low-mass star disks 
and 4/5 HAEBE star disks showing signs of crystalline 
silicates. Our sample is small, however, and may not be 
representative. Observations of a larger sample of low- 
mass stars are needed to quantify and understand the 
differences in the composition of dust around low- versus 
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Fig. 4. — SEDs and 10 /im spectra for disks around intcrmcdiatc-mass HAEBE stars where the silicate 10 /im band is in emission, 
otherwise the same as Figure 151 



intermediate-mass stars. 

5.3. Featureless Spectra 

The spectra and SEDs for debris disks around low- 
and intermediate-mass stars are shown in Figures and 
El respectively. No silicate emission was seen toward 
these disks — their spectra are predominantly photo- 
spheric, with a few exhibiting dust continua. The lack of 
silicate emission may indicate either a paucity of small 
grains (a < 10 /^m), or a dearth of optically thin ma- 



terial inside of ~10 AU (T > 200 K). The absence of a 
substantial continuum at 10 /im is consistent with the 
SEDs for HD 102647, HD 17925, HD 22049 and HD 
216803, which show very low IR and millimeter excesses 
(L/fl/L* < 10~ 3 ) with Xonset > 60 fim. This suggests 
that the debris disks have lost most of their (warm, 
small) circumstellar silicate grains. 49 Ceti, HR 4796A, 
and HD 233517 9 , however, still possess substantial far- 
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Fig. 5. — SEDs including 10 fj.m spectra for "debris" disks around 
low-mass stars. Dashed lines represent the stellar photosphere. 
The inset shows an enlargement of the 8-13 /im region with linear 
scaling on both axes. 



IR and millimeter excesses. This substantial amount 
of cold dust combined with the absence of silicate emis- 
sion suggest the possible existence of large gaps in these 
disks at radii close to the star, similar to those modeled 
in the debris disks a round e Eri, Vega, and HR 4796A 
ijOzernov et alJl2000t iWvatt et al.l[l999t lWvattjl2003|) . 

5.4. Evolution of Silicate Spectra 

The SEDs and silicate spectra of the sample indicate 
a possible evolutionary trend; silicate absorption is ob- 
served toward embedded protostars, strong silicate emis- 
sion arises from objects with optically thick protoplane- 
tary disks and nearly photospheric spectra are observed 
toward optically thin debris disks. The details of this 
trend are discussed below for both low- and intermediate- 

9 There is some evidence that HD 233517 m ay be a Li-rich giant 
and n ot a young star with a debris disk. Sec Balachandran et al. 
12000) for details. 
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Fig. 6. — SEDs including 10 fim spectra for "debris" disks around 
intermediate-mass stars, otherwise the same as Figure 



mass stars using a representative subset of our sample as 
shown in Figures [7| and |SJ 

Figure \7\ shows the evolutionary sequence for material 
around low-mass stars. The SED for IRAS 04239+2438 
indicates an embedded protostar, and the 8-13 /im spec- 
trum depicts deep silicate absorption (similar to 7 other 
embedded sources in our sample). The SED of LkCa 15 
is fit by an optically revealed star and near- to far-IR 
excess from an opt ically thick protoplanetary disk (c.f., 
IChiang et al.l 12001). and strong silicate emission is ob- 
served toward this source. Similar SEDs and spectra are 
seen for 3 other sources. The mid-IR spectrum of IRAS 
04295+2251 shows contributions from both silicate emis- 
sion, peaking near 10.5 /im, and silicate absorption near 
9.5 /zm and the SED suggests that this source (and 5 
others) may be intermediary between IRAS 04239+2438 
and LkCa 15, although the exact geometry is not yet 
well understood. As discussed above, the SED of Hen 3- 
600A shows no near-IR excess and very little far-IR flux 
and silicate emission peaking characteristic of forsterite 
(peaking near 11.3 /im), indicating substantial dust pro- 
cessing in this disk. Finally, HD 17925 is characteristic 
of 4 sources in our sample with SEDs consistent with a 
stellar photospheres (out to 60 /im) toward which no sil- 
icate feature is observed, indicating an absence of warm, 
small circumstellar grains, and longer wavelength excess. 

A similar evolutionary sequence is shown for the 
intermediate-mass stars (Figure [SJ, from silicate absorp- 
tion toward the fully embedded Becklin Neugebauer ob- 
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Fig. 7. — Evolution of the silicate spectrum for low-mass stars. The left panel displays the SEDs for selected sources in units of Log 
AF^, with the 10 fim spectra included for reference. The right panel shows an enlargement of the 8—13 /im region with linear scaling on 
both axes. Symbols and lines are as described in Figures I lltil The evolutionary sequence for low-mass stars is quite similar to that for 
intermediate-mass stars (Mccus & Waclkcns 1999). 



ject (BN) to silicate emission in the three optically thick 
disks and a nearly photospheric spectrum toward the 
optically thin disk encircling HD 102647. No sources 
with absorption+emission spectra were seen for the 
intermediate-mass stars. The sample, however, consisted 
of only three embedded protostars and there are large 
differences in morphology for these absorption spectra. 
Between the spectra of MWC 480, HD 163296 and HD 



179218, the silicate emission feature seems to evolve from 
amorphous olivine peaking near 9.5 /xm to crystalline 
forsterite, peaking at 11.3 /xm, suggesting dust process- 
ing, but there appears to be no corresponding evolution 
of the SEDs. Furthermore, MWC 480 shows silicate 
emission peaking near 10.5 /jm, indicative of crystalline 
cnstatite, not forsterite. 
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Fig. 8. — Evolution of the silicate spectrum for intermediate-mas s stars. Otherwise the sam e as Figure The evolutionary sequence is 
similar to that observed for intermediate-mass stars in ISO studies IMeeus fc WaelkenJn^99l) . 



6. QUANTIFYING SILICATE SPECTRA 

Having presented our spectra and discussed their over- 
all properties qualitatively, we now embark on a quan- 
titative analysis of the spectra. Detailed interpretation 
of silicate emission/ absorption features necessitates an 
understanding of the optical depth of the circumstel- 
lar envelope or the vertical structure of the circumstel- 
lar disk. Because our sample encompasses a variety of 
complex circumstellar environments ranging from deeply 



self-embedded young stars to optically thin debris disks, 
such detailed modeling will not be presented here. In- 
stead, we investigate empirical correlations between sev- 
eral properties of the silicate features and attempt to re- 
late the observed trends to circumstellar properties and 
dust size and composition. In §6.1, we begin by evaluat- 
ing the impact of the source morphology on the structure 
of the observed silicate feature by investigating correla- 
tions between the spectral shape and SED properties. 
Next, we explore the relationship between the FWHM 
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of the 8-13 fim emission/ absorption feature and the cen- 
tral wavelength of the peak/dip (§6.2), which are con- 
nected to optical depth and dust composition. The re- 
lationship between the 9.8 and 11.3 /zm fluxes and the 
peak-to-continuum ratio is indicative of grain size as well 
as crystallization fraction, as demonst rated by previous 
studies of HAEBE and T Tauri stars (jPrzygodda et all 
120031 Ivan Boekel et al.ll2003|) . Thus, relationships be- 
tween spectral shape and dust composition and grain size 
will be explored in §6.3-6.4. By applying these methods 
of analysis to our spectra, we hope to reveal more subtle 
evolutionary trends than are apparent from simple visual 
inspection of the data. 

6.1. Silicate feature strength and SED Properties 

Based on Figures and |H1 we have argued for evidence 
of an evolutionary sequence in the 10 /xm feature from 
absorption to mixed absorption and emission to pure 
emission, which may correspond to an evolutionary se- 
quence in SED morphology. To explore this relationship 
more quantitatively we investigate possible correlations 
between 10 /jm spectral properties and SED slope or IR 
e xcess luminosity. 

Lada & Wilking (1984) define a classification sys- 
tem using the SED slope, or 2-25 /im spectral index 
(o!2 -25 ^m)' The spectral index is negative for class I 
sources (embedded stars), positive in the range of 0-1 
for class II sources (optically thick disks), and positive 
and large (~3) for class III sources (optically thin disks). 
In Figure OH we explore possible correlations between the 
feature strength and the spectral index for spectra with 
pure emission or absorption features. Spectra showing 
evidence of mixed absorption and emission are not in- 
cluded. The feature strengths (F^ 1 ) used here are cal- 
culated from the normalized spectra (F= (F„ — F c )/F c ), 
by averaging over a region 0.2 (im wide centered on the 
largest absolute value (|F^ X |) in the normalized spectrum 
(see Table EJl ■ The spectra were smoothed before fitting 
to remove large fluctuations due to noise, which may have 
a large effect on locating the maximum. The quoted 
errors include the effects of spectral S/N (the standard 
deviation of the mean) as well as the effects of the contin- 
uum normalization method (the errors in determination 
of the endpoint fluxes). Error bars are particularly large 
for GM Aur, for which the spectrum possessed low S/N 
near the endpoints used for continuum normalization. 
For this reason, the error bars for GM Aur shown in all 
figures have been truncated at 1/7 of the full error. The 
correlation is evaluated using a x 2 linear least-squares 
regression method. During the calculation of the regres- 
sions, each data point is weighted by the inverse square 
of the error in x and y, so the effects of data with large 
errors are minimized. One should note that this method 
may reduce the weighting of T Tauri star disks, as their 
spectra for the most part have lower S/N than do spectra 
of disks around HAEBE stars. 

There is a strong inverse correlation between the 
strength of the absorption features and the spectral in- 
dex. As the spectral index increases, moving from class 
I toward the realm of class II sources, the strength of ab- 
sorption features becomes shallower. Furthermore, 4/5 
of the complex spectra (absorption+emission) , which are 
not included in Figure [51 possess spectral indices which 
are at the low end of the distribution (<0.5), suggesting 
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Fig. 9. — Correlations with the SED. This plot shows the cor- 
relation between the strength of the silicate feature (F\ 1 ; defined 
in Table 151 and the 2—25 /an spectral index (0:2— ihum)- The open 
symbols represent absorption spectra of solar-mass (squares) and 
intermediate-mass (triangles) stars in our sample. The solid sym- 
bols represent emission spectra. Spectra that can be fit with both 
emission and absorption are omitted. The error bars for Q2— 2Zum 
include the photometric errors as quoted in the literature. Fluxes 
are calculated by averaging the continuum normalized spectrum 
over 0.2 fim. The data are weighted by the inverse square of the 
errors shown, which include the effects of spectral S/N and the 
removal of the continuum. For clarity, the error bars for GM Aur 
have been truncated at 1/7 of the full error (see text). The strength 
of the absorption features appears to be strongly correlated with 
the spectral index (r = 0.75, 1.3% probability of being drawn from 
random distribution), as indicated by the dashed line. In contrast, 
the strength of the emission features is only weakly correlated with 
the spectral index, if at all; r = 0.29, with a 40% probability that 
the data are drawn from a random distribution. 



that they may be in an intermediate stage between class 
I and class II. This trend indicates that the strength of 
silicate absorption is primarily a function of the optical 
depth and not dust composition/size. In contrast, there 
is no compelling correlation between strength and spec- 
tral index for the emission sources in Figure 

Although the strength of the silicate emission features 
not related to the spectral index, it may still be de- 
pendent on or related to the disk morphology. In or- 
der to determine the relationship between the strength 
of silicate emission and disk morphology, the ratio of 
the IR to stellar luminosity was calculated. L/jj/L* 
is related to the optical depth of the disk material 
for optically thin emission. For optically thick emis- 
sion, however, hm/h* is indicative of the disk geome- 
try, tracing the fraction of light intercepted by the disk 
(Bac kman fc ParesceH993j) . In this case, disk flaring and 
orientation both affect the IR-to-stellar luminosity ratio, 
with a maximum value of L/ij/L* of 1/2 for a flared 
disk, in which the thickness increases with distance from 
the star (jKenvon fc Hartmannl lT987). and 1/4 for a flat 
optically thick disk. IR-to-stellar luminosity ratios for 
the disks in our sample, calculated as described in §4, 
are presented in Table [3] The optically thin debris disks 
possess L/fl/L* <0.1 and do not show evidence of silicate 
emission. We find generally higher values of L/^/L* for 
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the optically thick class II disks. These ratios are always 
greater than 0.1 with a maximum of ~0.5. The strengths 
of the emission features are not found to be correlated 
with L/fl/L*. As L/fl/L* is affected by both disk flaring 
and orientation, we can only say that the combined effect 
of the two is not clearly related to the strength of the 10 
/jm silicate emission feature. 

In summary, the strength of silicate absorption can be 
tied to the self-embeddedness of the star and is probing 
primarily the optical depth of the absorbing dust, the 
strength of the emission features appears to trace grain 
properties (e.g., grain growth and/or crystallinity), as 
opposed to optical depth, and these properties do not 
appear to be directly connected to evo lutionary stage, as 
traced by the lLada fc Wilkinel fl984jj) classification sys- 
tem, or to disk morphology. This is consistent with a 
scenario in which the observed emission arises from an 
optically thin region on disk surfaces. 

6.2. Width and Central Wavelength of the Silicate 
Feature 

For 6 lines of sight in the ISM, iBowev &: Adams"or] 
( 2002) find that the FWHM of the observed silicate ab- 
sorption decreases as the features shift to longer wave- 
length. They also find that these parameters are strongly 
related to the stellar environment. Three spectra probing 
the diffuse medium showed silicate absorption features 
with minima at ~9.8 fjxn. and FWHM of ~2.5 /im, typi- 
cal of amorphous olivine. The others, within the Taurus 
molecular cloud, possess silicate absorption features near 
9.6 /im with FWHM of 3.3 /im, typical of amorphous py- 
roxene. The authors suggest that this trend represents 
an evolution from olivine-dominated silicates in the dif- 
fuse ISM to pyroxene-dominated silicates in circumstellar 
environments. 

In order to search for trends in the silicate composition 
as a function of circumstellar environment, we perform 
similar calculations for the silicate features observed to- 
ward our sample of 17 embedded YSOs and 10 optically 
thick disks. Sources with photospheric spectra are not 
included. As discussed above, the strength (and shape) 
of silicate absorption features is a function of the op- 
tical depth. The central wavelength of the absorption, 
however, is dependent on the dust composition. In con- 
trast, the strength and width of silicate emission fea- 
tures (if optically thin) is almost solely dependent on 
the dust cross section, which is affected by both grain 
composition and size. The peak of the emission is also 
believ ed to be a function of b oth grain size and compo- 
sition ijvan Boekel et al.ll2003() . 

The silicate feature in each normalized spectrum (F= 
(F„-F c )/F c ) was fit by a gaussian to find the FWHM, 
central wavelength (A2), and the normalized flux (F>, 2 ) 
corresponding to that central wavelength. To account for 
the structured nature of the silicate features observed, 
we also record the wavelength (Ai) corresponding to the 
largest absolute value (|F^j |; used in §6.1) in the normal- 
ized spectrum. This provides a better characterization in 
the case of dust with large crystalline-to-amorphous sili- 
cate emission ratios (e.g., for HD 179218) where F\ 2 lies 
near the center of the feature at A2 ~10.5 /im, but F\ t 
successfully traces the stronger flux near the forsterite 
peak (Ai ~11 /im). The results of these calculations are 
shown in Tabic |5| 
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Silicate spectral parameters 
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-0.15 
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9.61 


10.02 


1.33 


1.34 


2.48 


cm 


HD 179218 


11.01 


10.48 


0.56 


0.53 


2.39 


cm 


WW Vul 


9.84 


9.92 


1.65 


1.63 


2.33 


cm 


HD 184761 


10.06 


10.07 


1.45 


1.46 


2.20 


cm 



Note. — Ai and Fi indicate central wavelengths and feature 
strengths calculated by fitting a gaussian to the normalized spectrum. 
A2 and F2 indicate central wavelengths and feature strengths calcu- 
lated by finding the largest absolute value of the normalized flux. In 
all cases, FWHM are found from gaussian fits. 



a The observed spectrum can be fit by absorption + emission. 



This analysis is not as straight-forward for the spec- 
tra that show evidence of both emission and absorp- 
tion. If we fit the absorption and emission components 
individually for the complex spectra, as shown in Ta- 
ble 03 we find small FWHM for the individual features, 
which are shifted to shorter/longer wavelengths for ab- 
sorption/emission. For all 5 of the complex spectra we fit 
absorption components with Ai,2 < 10 /an and emission 
components with Ai,2 > 10 /im. Although these complex 
spectra represent an interesting new class of sources, the 
emission and absorption features clearly can not be de- 
convolved empirically. Source geometry must be taken 
into account in the interpretation of these spectra. For 
these reasons, we include only pure absorption and emis- 
sion features in the rest of our analyses. 

In Figure ITU1 the central wavelengths (Ai and A2) of 
the emission/absorption maxima (as described by F^ x 
and F^ 2 ) are plotted with respect to the FWHM. Al- 
though a subset of the pure absorption or emission fea- 
tures are clustered about 9.75±0.15 /Ltm, the sample en- 
compasses a large range of central wavelength (~9-ll 
/im) and FWHM (~ 1.7-3 /im), suggesting a range of 
source environments. There appears to be no strong cor- 
relation between A 1.2 and FWHM in the sample as a 
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Fig. 10. — Plot of the wavelength of the emission/absorption 
maxima (A1.2) and the width of the feature (FWHM), as presented 
m Table y Ai and A 2 are represented by connected pairs. Sym- 
bols are as described in Figure [5] Neither the sample of silicate 
absorption spectra nor that of silicate emission spectra show sig- 
nificant correlation between Ai,2 and FWHM. The two samples do 
show different trends, however, with the emission spectra spanning 
a larger range in central wavelength and the absorption spectra 
spanning a larger range in FWHM. 



whole [r = 0.04), nor in either sub-population of silicate 
absorption and emission spectra. The two samples do 
show different trends, however, with the emission spec- 
tra spanning a larger range of central wavelength and the 
absorption spectra spanning a larger range of FWHM. To 
address the nature of these trends in more detail, we now 
discuss the absorption and emission subsamples in turn. 

The silicate absorption features show little variation 
in central wavelength (9.6±0.2 /im), which is consistent 
wi th amorphous olivines and pyroxenes. The trend seen 
by Bowcv & Adamson ( 2002) is not seen, suggesting that 
there is no clear correlation within this sample of class I 
objects between the grain composition (traced by 
and optical depth (traced by FWHM). The FWHM of the 
absorption features varies significantly (~1.8-3.2 /<m), 
indicating a range of optical depths, indicative of varying 
envelope thickness or morphology. There are three out- 
lying absorption sources possessing large FWHM (IRAS 
04287+1801, IRAS 04381+2540, and Mon R2 IRS3). 
These were previously identified in §5.1 as possessing ad- 
ditional shallow feature near 11.2 /im and the very large 
FWHM may be merely a result of the blending of the 
11.2 jLim ice feature with the silicate feature. 

The silicate emission features in general peak at longer 
central wavelength (10.2+0.3 /im) in comparison to the 
embedded sources, indicating either grain growth or 
modified composition. The silicate emission features also 
show larger variations in the central wavelength, but lit- 
tle variation in the FWHM. The FWHM for the emis- 
sion features is largely dependent on the grain cross sec- 
tion and, thus, on both the grain size and composition. 
The range in FWHM of roughly 2 to 2.5 /im is consis- 
tent with primarily amorphou s silicate grains of 0.1-2 /im 
sizes i|van Boekel et al.ll2003j) . The variations in FWHM 



are also consistent with those expected from the addition 
of crystalline silicates. Indeed, the 2 emission spectra 
with the largest peak wavelengths (~ 11.2 /im) are those 
previously identified as possessing strong crystalline sili- 
cate emission, Hen 3-600A and HD 179218. The emission 
spectrum from MWC 480 is also an outlier, peaking at 
10.3-10.5 /xm, which is consistent with crystalline en- 
statitc. 

Thus, sources with silicate absorption features have 
central wavelengths and FWHM generally similar to 
amorphous olivine or pyroxene, while those with silicate 
emission show increased peak wavelengths and variations 
in FWHM characteristic of more processed dust. The re- 
lationship between the silicate emission feature and size 
and compositional variations of the dust will be explored 
in more detail below. We will not discuss the absorption 
or the absorption-plus-emission sources any further. 

6.3. Spectral Color-Color Diagrams: Grain Composition 

Now we attempt to use fluxes derived from the con- 
tinuum normalized silicate emission features as probes 
of grain composition, following the methods used fo r 
ISO spectra of HAEBE stars l|Bouwman et alJ f20011. 
In particular, we will compare fluxes at the nominal 
wavelengths of emission from amorphous and crystalline 
olivine, 9.8 and 11.3 /im, and ratios of 9.8-to-8.6 versus 
9.8-to-11.3 /j,m fluxes in an attempt to assess the probe 
the crystallization mechanism. Figure ^] displays flux 
and color-color plots for the silicate emission features ob- 
served toward our sample. The fluxes are calculated by 
averaging over 0.2 /im. The errors are calculated as de- 
scribed in §6.1 and include the affects of S/N and the 
continuum subtraction method. 

In Figure ITTk. the relationship between the blue and 
red slopes is examined. Assuming the emission arises pri- 
marily from amorphous olivine (peaking near 9.8 /im), 
the slopes on the long and short wavelength side of the 
normalized spectra can be approximated by F9.8/F11.3 
and F9.8/F8.6, respectively. The dashed line indicates the 
correlation (r = 0.68) between the blue and red slopes 
for the observed silicate emission features. This trend im- 
plies that decreasing Fg.g/Fg.g (blue slopes) correspond 
to decreasing F9.8/F11.3 (red slopes), indicating that the 
entire feature becomes flatter as emission at 11.3 /im be- 
comes more prominent. The two sources with silicate 
features that peak above 11 /im, HD 179218 and Hen 3- 
600A, also fall along this trend, occupying the lower left 
corner of the plot. Although the variations in the struc- 
ture of the silicate emission has often been attributed 
to variations in crystalline fraction, this trend is consis- 
tent with grain compositions dominated by amorphous 
olivine, with variations due to changes in grain size (e.g., 
Ivan Boekel eTaLll200^ . 

If the above trend can be attributed to the conversion 
from dust dominated by amorphous olivine to that domi- 
nated by crystalline forsterite, then variations in the crys- 
tallinity of the dust in these sources should be evident in 
the relationship between the fluxes at these wavelengths. 
The fluxes corresponding to the peaks of the amorphous 
olivine (Fg.g) and crystalline forsterite features (Fn.3) for 
the normalized spectra are shown in Figure ITTb . There 
is a linear correlation between the fluxes at 11.3 and 9.8 
/im indicating that the flux at 11.3 /im does not increase 
at the expense of the flux at 9.8 /im, as we would expect 
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Fig. 11. — IR color-color diagrams for silicate emission, (a) 
The relationship between the slopes on the long (Fg.g/Fn.3) and 
short wavelength (Fg.g/Fg.6) side of the normalized spectra. A 
linear least-squares fit to the data is depicted with a dashed line, 
(Fg.s/Fg.s) = (0.7±0.2)x(F 9 . 8 /Fii. 3 ) + (0.6±0.2), and character- 
ized by a correlation coefficient of r = 0.67. (b) Fluxes correspond- 
ing to amorphous olivine (Fg.g) and crystalline forsterite (F11.3) 
in the normalized spectra. The linear fit to these data (dashed 
line, r = 0.85), Fg.g = (3.0±0.6) X F11.3 - (3.0±1.0), indicates 
that F11.3 does not increase at the expense of Fg.g. Fluxes and 
corresponding error bars are calculated as in Figure 151 



from crystallization of amorphous olivines (emitting near 
9.8 /im). Additionally, the spectra with dominant emis- 
sion at 11.3 /im, are weaker overall and fall in the lower 
left corner of the plot. This suggests that the crystalline 
silicates observed toward these sources either are not pro- 
duced from the same reservoir of small amorphous grains 
that emits at 9.8 /im, or that variation in the shape of 
the silicate feature is not dominated by crystallization. 

The above analysis demonstrates that the 9.8-to-8.6 
and 9.8-to-11.3 colors are correlated within the observed 
silicate emission features. The features become flatter as 
the relative emission at 11.3 /im increases, which is con- 
sistent with increased grain sizes. Furthermore, the flux 
at 11.3 /im does not decrease at the expense of that at 



9.8 /mi, indicating that if the flattening of the feature is 
due to crystallization, then the crystalline forsterite emit- 
ting at 11.3 /im are not crystallized from the amorphous 
silicate population emitting at 9.8 /im. 

6.4. Strength and Shape: Grain Size 

Here we discuss the possibility that the changes in the 
silicate emission feature shape discussed above are due 
to changes in grain size, rather than crystallinity. Al- 
though observations of broadened 8-13 /im silicate fea- 
tures with emission observed at 11.3 /zm were long be- 
lieved to be an clear indication of the presence of crys- 
talline fors terite, the analysis prese nted in §6.3 and recent 
studies by ( van Bockel et al. 2003) suggest that such fea- 
tures are actually representative of increased grain size. 
These studies show that the mid-IR absorption coeffi- 
cient of ~0.1 /im-sized amorphous olivine grains pos- 
sesses a triangular shape similar to that observed in the 
diffuse ISM. Due to a general decrease in the absorption 
cross section, a population of larger amorphous olivine 
grains (~2.0 /zm) has an absorption coefficient that is 
weaker by a factor of ^2.3 and flattened between 9.5 
and 12 /im. Examination of the emission spectra of 
12 HAEBE stars showed that there is a correlation be- 
tween the feature strength (peak flux of the normalized 
spectrum) and shape (as indicated by the ratio of the 
flux at 11.2 to 8.9 yttm), with broader emission features 
being weaker as shown in Figure IT2*fa (open triangles). 
iPrzveodda et alJ l)2003|) found a similar trend for a sam- 
ple of disks surrounding low-mass stars (open squares). 

We perform a similar analysis for our sample of T Tauri 
and HAEBE stars, which are shown as solid squares and 
triangles in Figure 112b . Fluxes and errors are calcu- 
lated as described in §6.3. The entire sample of both 
T Tauri and HAEBE star spectra can be fit by a single 
trend. This trend is consistent with variations in grain 
size among the sample, but we can not rule out the pos- 
sibility that the changes in feature shape and strength 
result from an increased crystallization fraction (§6.3). 
The two sources with spectra consistent with forsterite 
emission at 11.3 /im do not deviate significantly from this 
trend. It is therefore reasonable to think of the strong, 
narrow and the broad, weak silicate emission features 
as arising from "unprocessed" and "processed" dust, re- 
spectively. The contributions from dust size and com- 
position can thus only be separated when emission from 
crystalline silicates are clearly seen in the spectrum as 
an isolated feature at 11.3 /im. 

Although generally used as a measure of the crystal- 
lization fraction in the emitting silicate dust (e.g., the 
"9.8-to-11.3 phase transition"), the wavelength of maxi- 
mum emission may also be affected by grain size. In Fig- 
ure the maximum emission in the normalized spec- 
trum Fx! is plotted with respect to the wavelength cor- 
responding to the maximum Ai. Figure indicates 
that for our sample as the feature maximum shifts to 
longer wavelengths, the strength of the emission feature 
decreases. This trend is consistent with the effects of 
grain growth. As the grain size increases, the feature 
flattens and the maximum decreases in flux and shifts 
to longer wavelength. Although, the two sources with 
strong crystalline forsterite features (HD 179218 and Hen 
3-600A) have weaker peak flux and show up at the bot- 
tom right in Figure 112b . in general moving along the 



20 



Kessler-Silacci et al. 



CO 
02 




1.0 1.5 2.0 2.5 3.0 3.5 



A, 



-< 



2.5 
2.0 

1 .5 

1 .0 

0.5 
0.0 




9.0 9.5 10.0 10.5 1 1 .0 1 1 .5 



Fig. 12. — Emission feature characteristics. In both plots, the 
solid squares and triangles represent the observations of T Tauri 
and HAEBE stars in our sample, (a) The correlation between fea- 
ture shape (F11.3/F9.8) and strength (F\ 1 ), as noted previously for 
HAEBE (van Bockcl ct al . 2003, open triangles) and T Tauri stars 
IPrzygodda et al.M20oT open squares). The dashed line is a linear 
fit to the entire sample, Fii. 3 /Fg.8=(1.41±0.07)xF Al -(0.27±0.03), 
with a correlation coefficient of r = 0.82. (b) The strength of the 
silicate feature decreases with decreasing peak wavelength (Ai) as, 
F Al = (-1.0±0.1) X Ai + (11±1), with r = 0.78. Both trends are 
consistent with grain growth (see text). Fluxes and corresponding 
error bars are calculated as in Figure l9l 



trend from the upper left to lower right of Figure 112b 
does not coincide with the appearance of emission from 
crystalline forsterite. HD 163296, which has a visible 
bump near 11.3 /im, has a silicate feature which is still 
quite strong (in the upper left of the plot). Further- 
more, the broad silicate emission feature observed toward 
AA Tau, which shows no evidence of emission from crys- 
talline forsterite, is quite weak and lies in the middle of 
the trend. Therefore, the shift in the wavelength of max- 
imum emission in the 10 /im region, or the "9.8-to-11.3 
fim phase transition" , appears to be largely dependent 
grain size variations and thus is in fact not a good indi- 



cator of crystallization fraction. 

If the trends depicted in Figure El primarily indicate a 
change in the grain size, then they suggest a progressive 
removal of small grains (< 2 /im) from the warm (T rj 
300K) disk surface probed by the 10 /im feature. Dust 
settling would preferentially remove larger grains from 
the disk surface. The absence of small grains could be 
related to the removal of dust through gap clearing, but if 
so we would expect that GM Aur, which shows evidence 
of a substantial gap, to display a broadened feature in- 
dicative of larger grains. The flat, broadened spectra in 
the bottom right of Figure 112b are therefore most likely 
a result of grain growth in the surfaces of these disks. 

7. CONCLUSIONS AND FURTHER WORK 

8-13 fim spectra were obtained for 34 low- and 
intermediate-mass YSOs. The evolutionary sequence 
seen i n ISO observations of high/intcrmediate-mass 
YSOs <E 

eeus fc Waelkenall999|) . was confirmed and ex- 
tended here to ~solar-mass stars. Amorphous silicate 
absorption is observed toward young, embedded sources. 
In 5 of 17 embedded sources, complex spectra that may 
be interpreted as a combination of silicate absorption and 
emission are observed. These spectra may represent an 
interesting new class of sources that are in between Class 
I and Class II, but require physical models to decon- 
volve the emission and absorption. Silicates are observed 
in emission toward optically revealed stars. The shape 
of these emission features varies from strong triangular 
emission peaked near 9.8 /im, to flat, broad emission with 
a central wavelength between 10 and 11 /im. Finally, 
silicate emission was absent for the debris disks in our 
sample. 

To aid in the interpretation of the silicate features, we 
explore the connections between silicate feature proper- 
ties with stellar evolution through comparison with the 
spectral energy distributions. We find that the strength 
of 10 /im silicate absorption is strongly correlated with 
the 2-25 /im spectral index. Class I sources with complex 
absorption+emission features have small spectral indices 
when compared to those with pure absorption, suggest- 
ing that the former are less embedded. These results 
indicate that the silicate absorption feature is indeed a 
good probe of the circumstellar environment, tracing the 
self-embeddedness of the star. In contrast, the strength 
of the silicate emission feature does not appear to be cor- 
related with the 2-25 /zm spectral index (evolutionary 
stage) or the IR-to-stellar luminosity ratio (disk geome- 
try) and thus probes grain properties rather than source 
morphology. 

Although the strength of the silicate absorption fea- 
tures is dependent on the circumstellar environment, the 
central wavelength of the emission should be sensitive 
to the grain composition. We therefore attempted to 
quantify the effects of grain properties, such as size and 
composition, on the shape and/or central wavelength of 
all of the silicate features observed in our sample. In 
general, the features are clustered around 9.7 /im, corre- 
sponding to amorphous olivine, with FWHM near ~2.3 
/im. A direct correlation between FWHM and central 
wavelength is seen for neither silicate emission nor ab- 
sorption features. The pure silicate absorption features 
exhibit a large range of FWHM, likely related to differ- 
ences in optical depth, and narrow range of central wave- 
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lengths, indicating limited variation in grain size and/or 
composition from that of the ISM. A larger scatter and 
higher mean central wavelength is found for the emission 
spectra, likely related to compositional variations within 
the sample (i.e., fraction of amorphous vs. crystalline 
silicates), indicating increased dust processing. 

As the emission features are more sensitive to grain 
properties, we further analyzed the structure of the emis- 
sion features, finding: 

- The shape of the silicate emission of most T Tauri 
stars was similar to that previously observed in the 
diffuse ISM, in contrast to the structured emission 
feature observed for most HAEBE stars. This sug- 
gests an absence of small (< 10 /J.m), warm (T«300 
K) crystalline silicates in these T Tauri star disks, 
but colder or larger crystalline silicate grains may 
still be present and, if so, will be visible through far- 
infrared lattice modes accessible from space born 
observations. 

- A correlation between the peak shape as traced 
by F11.3/F9.8 and feature strength for 5 T Tauri 
and 5 HAEBE in our sample was also found. This 
is con sistent with that seen for several HAEBE 
stars frvan Boekel et al.l l2003|) and T Tauri stars 
l|Przvgodda et al]l2003j) . which are attributed to 
changes in grain size (from 0.1-2 /jm). The entire 
sample can be fit best with a single trend, includ- 
ing spectra with strong indications of crystalline 
forsterite. The wavelength of maximum emission 
was also found to be correlated with peak strength 
for our sample. This indicates that interpreta- 
tion of the "9.8-to-11.3 /jm phase transition" must 
therefore include analysis of grain size in addition 
to crystallization fraction. 

- We also demonstrate that the flux at 11.3 fim 
does not increase at the expense of that at 9.8 
fim. If crystallization is dominating the transition, 
then this indicates that the crystalline forsterite 
emitting near 11.3 /im was likely not formed from 
the population of amorphous olivine emitting near 
9.8 /xm, suggesting a more complex crystallization 
mechanism. Another interpretation is that the re- 
lationship between the 11.3 and 9.8 /im fluxes is 
dominated by grain size effects and not crystalliza- 
tion. 

10 Details of the C2D and FEPS programs ca n be 
found at |http://w ww.peggysue,as. utexas.edu/SIRTF and 



Silicate emission features in solar-mass T Tauri stars 
have primarily been observed in the 10 /jm atmospheric 
window (as in this study). Due to sensitivity limitations, 
ISO studies of silicates focused on intermediate or high- 
mass stars. The Spitzer Space Telescope will be used 
to expand these studies for a large sample of low-mass 
stars, ranging from embedded protostars to optically thin 
disks, via the combination of various GTO programs and 
the C2D (T dust > 1) and FEPS (T dust < 1) Legacy pro- 
grams 10 . In this manner, the study of the crystalline 
and amorphous silicates in disks, will be expanded to in- 
clude disks around low-mass, sun-like stars, creating a 
database analogous to ISO studies of high/intermediate- 
mass stars. These studies will include longer wavelength 
silicate bands, which depend intricately on the coordina- 
tion of the silicon atoms and can provide more specific 
information about the minerals present, including the 
Mg/Fe ratio. The mineral content of the grains in each 
stage of evolution can thus be directly compared to that 
of meteorites, asteroids and comets. 
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APPENDIX 

VARIABILITY OF 8-13 fiU SILICATE EMISSION SPECTRA 

Long term monitoring of several sources, including DG Tau l|Wooden et al.ll200fl and HD 45677 (|Sitko et alJ ll994). 
has revealed dramatic variations in the strength and shape of silicate features over periods of months to years. Ob- 
servations of DG Tau revealed that although the dust continuum remained fairly constant, the silicate feature evolved 
from nearly featureless continuum to emission from crystalline and amorphous silicates to e mission from primarily 
amorphous silicates and, finally, to absorption from amorphous silicates between 1996 and 2000. lWooden~ al. ( 2000) 
indicate that these variations in the silicate feature are due to dust in the op tically thin disk surface, and that the 
continuum is produced by the underlying disk midplane (e.g.. lCalvet et af]ll992j) . For HD 45677, the flux of the silicate 
emission featu re was observed to decrease by —30% between 1980 and 1992, coincident with an increase in the stellar 
flux of -60% ( Sitko et al. 1994). The authors suggest a scenario in which long-term (decade-long) variations in UV 
and optical flux are due to a clear hole moving into the line of sight, which is possibly coupled with a loss of small 
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Fig. A13. — Searches for variability. The left panel shows the silicate emission toward HD 163296 observed on 1999 August 23 (solid line) 
and 2000 June 20 (dashed line). The right panel shows the silicate emission observed toward LkCa 15 on 1999 November 30 (solid line) and 
2000 November 8 (dashed line). The spectra from two observing runs are overlaid, showing that there is no significant variability (above 
the error) in the shape of the silicate emission feature for these two sources on a yearly timescale. The continuum was not subtracted for 
these spectra. 



grains, as indicated by the change in flux of the silicate feature. 

Multiple spectra were obtained toward a fe w sou rces in our sample, permitting study of the variability in shape 
and strength of the silicate emission. Figure IA13I shows the silicate emission toward HD 163296 and LkCa 15 for 
two observations taken approximately one year apart. Spectra from the two observing runs are overlaid, showing 
that they are consistent (within the errors). The doub le peaked structure se en in the spectrum of HD 163296 has 
been previou sly observed in s pectra obtained by ISO ijBouwman et al.1 120011 1996 October 14) and ground-based 
observations iSitko et alJll9 99, 1996 October 10) and these spectra do not differ significantly from our observations. 

HD 179218, Hen 3-600A and MWC 480 h ave also been observed on more than one occasion reported in the literature. 
The previous observations of HD 179218 l|Bouwman et all 1200 li ISO LWS, 1996 October 5) are virtually identical 
to our observations (2000 February 21). The spectra depict ve ry similar peak strengths and shapes. MWC 480, 
however, exhibits very different emission in the ISitko etldl ©99) (IRTF BAAS, 1996 October 14) study than m our 
spectrum (1999 November 30). In particular, the 11.3 /im feature (crystalline forsterite) is much less prominent in 
our observations and there is an additional peak at ~10.5 /im (near the position of crystalline enstatite). However, 
the (peak-to-contin uum) strength of the feature is the same to within ~15%. The observations of Hen 3-600A by 
Hond a et al.l {2003) (Subaru COMICS, 2001 December 27) also suggest strong evidence of variability when compared 
to our data (Figure 2000 February 21). Our spectrum shows much less emission at 9.2 /im and possibly 10.3 /jm 
(both attributed to enstatite), with the peak at 11.3 /im (forsterite) dominating the spectrum. There also appears to 
be an accompanying decrease in the continuum or in the strength of emission near 12.5 /im from 2000 to 2001. 

There is no obvious trend in the variability with respect to spectral parameters or silicate emission feature shape from 
these observations. HD 179218 and Hen 3-600A have similar emission features, ages and are in similar evolutionary 
stages (although their mass differs by a factor of 2). Yet the spectrum of Hen 3-600A is highly variable on 1 yr timcscalcs 
while the spectr um of HD 179218 has not changed in > 3 yr. Although HD 163296 and MWC 480 possessed similar 
spectra in 1996 l)Sitko et alJl!999» . the spectra for these sources are quite different from each other in 1999 (our data). 
MWC 480 has undergone changes in the spectrum that indicate conversion between the crystalline forms of forsterite 
(Mg2Si04) and enstatite (MgSiOs). This conversion, also seen in the H en 3-600A spectrum, is the opposite of the 
expected transition with time during the standard annealing process (see lGaillll998j) . 

Studies of the variability of silicate emission may provide valuable information about the mechanisms of grain 
growth and/or crystallization that affect the feature shape. However, variability of the silicate emission feature has 
been observed toward only a few objects, and rarely in a comprehensive program designed to monitor the accompanying 
variability of the star and/or SED. Dedicated long term monitoring of the silicate feature toward a larger sample of 
objects is required to understand the cause of silicate feature variability in disks. 

SED REFERENCES 

The SEDs presented in Figures I1I6I were compiled from photometric observations in the literature. The references 
for each source are listed below. 

49 Ce ti: T he SED inc ludes optical photometry fr om lHa ves fc L atham! (I1 975J). UBV pho t ometr y from iRvbkal 
<|1969ft and lStovl ill 9691) JHKL M photometry from lEiroa et all i200llTandlSvlvester et all l|1996l) . mid-IR pho- 
tometry from lThi et alJ 1)200 lj) , and IRAS observations from lBeichman et al.l l)1988|) . 



HD 17925: The SED includes optical photometry from lHaves fc Latham! l|1975f) . UBVRIJHKL photometry from 
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Johnson et all 11966ft . ICutispoto et all <|2001|) . lBessell 11990ft andlJofmson et etl.1 (|1968T> . N-b and photomet r y from 
Metchev et all fl2*003ft . and IRAS observations from IBeichman et alJ l)1988ft as presented in lHabing et alJ |2001). 



HD 22049: The SED includ es optica l phot ometry from lHaves fc Latham! 11975ft . Geneva 7-color ph^tome- 
try (U ,B1, B,B2,V1,V,G ) from iRufeneil lfl976ft . UBVRIJHKL photomet ry fr om I Johnson et all iTL966|) . iMendozal 
( 1970|) and lGlassI 11975ft. IRAS photome try from lBeichman et all l)1988ft and lHabing et al.l<j2001ft . and 800-1300 
/xm photometry from lChini et all l|199()). 

IRAS 04016+2610: The SED includes optical, J HKL , (sub) mil li meter photometry from iKenvon et all lTl 993) 
and re ferences therein, 1-2 ^m photometry from IPadgett et alJ 11999ft. mi d-IR photometry from IMvers eTlil] 
l|1987ft and this paper, submi llimeter photometry fromlYoung et al.H2003ft andlHogerheiide fc Sandelll (120001) . an d 
millimet er photome t ry fro m Moriartv-SchicvenetalJ 11994ft. iMotte fc Andrei 1)200 lft . iHogerheiide et alJ (|1997ft . 
lOhashi et ail l|1996ft . ISaito et all l|2001ft and lLucas et all l|2000ft . 

IRAS 04108+2803: The SED includes optical, J HKL, (sub) m i llimet er photometry from IKenvon et all (|1993fl 
and references therein , mid-IR photometry from iMvers et ; alj jl9$7ft and t his paper, IRAS photometry from 
I Beichman et alJ 11988ft. submill imeter photometry fro m lYomiget^Ll (12003^) . a nd millimeter photometry from 
iMoriartv-Schieven et all l)1994ft . IMotte fc Andrei l)2001ft and lOhashi et all l|1996^ . 



IRAS 04169+2702: The SED i ncludes K-band pho t ometr y from lTamura et all (11991ft. K-band, mi l limete r and 
submillimeter photometry fr omlBar sonv fc Kenyonl ((1992) , IRAS photometry f rom IB eichman jivLaU (tLQ SSft. sub - 
millimeter photometry from [Young et all (J2003), mi llimeter ph otometry from IMoriartv-Schieven et al.l (J1994) , 
IMotte fc Andrei (|2001^ . and radio photometry from ISaito et alJ l)2001ft . 

IRAS 04181+2654 A/B: The SEDs inclu de JHK band p h otome try from IKenvon et all l)1993ft and ref- 
erences therein, IRAS phot ometry from IBeichman et alJ l)1988ft . and submillimeter photometry from 
IMoriartv-Schieven et all l|l994ft . 

IRAS 04239+2436: The SED includes JHKL band p hotometry from IKenvon et all (11993ft and IPadgett et all 
(199 9[), IRAS photometry fromlBeichman et alJ 0.988), submil limeter photom e try fr o mlY oung et alJ (12003ft and 
Moriartv-Schieven et all l)1994ft . millimeter ph otometry from IMotte fc Andre! ()2001ft . lOhashi et alTl)1996ft . and 



s1uTcT^tHrtT((2"ob"l() . and radio photometry from lLucas et alJ (|200 



IRAS 04248+2612: The SED includes JHKL band photometry fr om IKenvon et al l (11993ft and IPadgett et alJ 
l|1999fl . N-ba nd photometry from this paper, IRAS photometry fr om |Beichma^^L^L|l| l 988ft . submillimeter pho- 
tome trv fromlYou ng et ah| 12003ft . IMoriartv-Schiev en et alJ l)1994ft . and lDent et alJ 11998ft. millimeter p hotometry 
from lMotte fc Andrei l)2001ft and lOhashi et alJ l)1996ft . and radio photometry from lLucas et all l)2000ft . 



IRAS 04264+2433: The SED includes JHKL band photometr y fromlKenvon et all 11993ft and references therein, 
N-ba nd photometry from this paper, IRAS photometry fro mlBe ichman et^dTi^988ft . submillimeter photometry 
from lYoung e t al. ( 2003), an d millimeter photometry from lMotte fc Andrei l)2001ft . 

Haro 6- 10 A/B: The SEDs inclu de BVRIJH and mid- IR photometry from IMvers et alJ 11987ft KLM pho- 



riaro b- lU A/ a: ine si^Ds includ e avtiUtl and mid- lit pnotometry trom IMvers et alJ li lVJSn. KLM prio- 
tometrv IKenvon et alJ l|1993ft an dlLeinert fc Haa^ dl9 89), IR AS photomet r y from IBeichman et alJ 1^)88), and 
(sub)millimeter photometry from lChandler et alJ 1 1998ft and IKenvon et alJ (l 993ft . 



IRAS 04287+1801: The SED includes JHKL M band photomet r y from IKenvon et alJ l)1993ft . N-band photom- 
et ry from this pap e r, IR A S photometry from IBe ichman et alJ iu988ft. and (sub)m i llimeter photomet r y from 

E handler fc Richerl (|2000ft. iMoriarty-Schieven et alJ l|1994ft . IMotte fc Andrei pOOlft . IHogerheiide et all (j 1 997ft . 
hashi et al.l l)199a) and lSaito et alJ 



IRAS 04295+2251: The SED includes JHKL and mid-IR photometry fr om IKenvon et alJ 11993ft and refer- 



ences therein, N-ban d photometry from this paper, IRAS photometry f rom |Belcliman^t^iT " () 1988ft . submillime- 
ter photometry from | Young e^l] ([2QQ3J and [ Moriartv-Schieven et alJ l|1994ll . and millimeter photometry from 
IMotte fc Arldrl (|2001ft and lOhashi et all (|1996ft . 

AA Tau: The S EP in cludes UBV photometry fromlJovl 11949ft mid-IR , photometry from lThi et all (|2001ft and 
| Metchev et al.l 12003ft. IRAS photometry from IBei chman et alJ 11988ft . and ( sub ) mill imet er photometry from 
Dutrev et alJ l| 1 996ft iBeckwith fc Sargent! ()1991ft . lBeckwith et alJ (|l99q) . and|^3 l|200l|) . 

LkCa 15: The SED inclu des UBVRIJHKL N photometry from lKenvon fc H artmann|l[l995) and references therein 



mid-IR photo metry from |Thi_el_al (2QQ1), IRAS photometry f rom iBgichman et alJ l)1988ft . and millimeter pho- 
tometry from lOsterloh fc Beckwithl l|1995ft . iDuvert et alJ l|2000ft . andlQil l|2001ft . 



IRAS 04381+2450: The SED inclu des JHKL and mid-IR photometry from IKenvon et all 11993ft and refer- 
ences t herein, IRAS photometry from | Bei chman et alJ119 88f). and (sub )mill imeter photometry fromrYoung et all 
l|2003l) . IMoriartv-Schieven et alJ l|T994ft . IHogerheiide fc SandeiH (|2000ft . and lHogerheiide et all l|1997ft . 
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IRAS 04489+3042: The SED includes JHKL and mid-IR photometry from iKenvon et alJ (119931) and references 
therein, N-band photometry from this paper, and millimeter photometry from iMot^^^Andrel 1)200 ID . 

GM Aur: The SED i nclud es UBV photometry from iBastian fc Mundtl l|1979|) . UBVRIJHK p hotometry from 
IKen yon fc Hartmannl l(1995|) . mid-IR photome try from |Th^O0l20^)7^I RAS photometry from |Bejchmane^l ' 
(1988|), and (sub)mill imeter photometry from lBeckwith fc Sargentl 1)1991(1 . iDutrev et all l)1996|) . iBeckwith et al 
and S (200 



MWC 480: The SED includ es optical p hotometry from lEimontas fc Sudziusl dT998^ and lHaunt fc Schrolll 
(I1974D. HKL photomet ry from lAllenl l)1973l) . mid-IR pho tometry fromlThi et all 120011) . IRAS phot ometry from 
Bcichma n et alJ (|1988|) . and millimeter photometry from (Mannings & Sargentl l|1997|) and Qi (200JJ). 

BN : The SED i ncludes near-IR pho tometry from iHillenbrand et all (I2001D andlDougados et all 1(19931) . mid-IR 
phot ometry fromlGezari et al.l lll998[K millimeter p hotometry from lPlambeck et al.l l|1995|) . and radio photometry 
from lFelli et ail II1993D and lMenten fc R.eidl 111993) . 

NGC 2024 I RS2: The SEP inclu des optical photometry from Hillen brand, unpublish e d. CCD, JHKL pho - 
tometry from | H aisch et aT| J2001f) . mid- and far-IR phot ometry from 



_ , - Grasdalcn (1974). Haisc h et alJ i2001). 

iThronson et alJ <|1978lll984t). millimeter photometry fromlWilson et al (119951) andlEisner fc Carpenter! <|2003|) . 
and radio photometry from lGaume et all <|1992|) . iRodri'guez et alJ lj2003f) . and lSnell fc Ballvl ljl986f) . 

Mon R2 IRS3: The S EP includes JHKLM photo metry fromlPreibisch et all ll2002T).lCarpenter et all l|1997|) . and 
iKoresko et all II1993D mid-I R photometry from iWalsh et alJ ll2001l) andlPersi et all (119961) . IRAS photometr y 
from iBeichman et alJ 119881), (sub)millimeter p hotometry from Uenness et alJ <|1995|T and lHenning et alJ ([l992), 
and radio photometry from lTofani et all l|1995|) . 

HP 2 33517: The SED includes opt ical photometry from ISvlveste r et alJ II1996D a nd iMiroshnichenko et alJ 
(I1996D. JHKLL'MN ph otometry from ISkinner et all lll995D andlSvlvester et alJ (I1996D. IRAS photometr y from 
IBeichman et al.l (|1988l) , and millimeter photometry from lPuncan et aTll|T990fr and ISvlvester et al.l l|2001|) . 



Hen 3-600A: The SEP includes B VRI photometry f rom iGregorio-Hetem et all (119921) and iTorres et alJ 
ll200l. .IHKLMN photom etry from IGeoffrav fc Moninl J2001D . IWebb et alJ (I1999T). IWeintraub et all 112000]) 
Uavawardhana et "all l|1999|) . and lMetchev et alJ <|2003|) .and IRAS photometry from lBeichman et alJ (|1988^ 



HP 102647: The SEP inc ludes op tical photometry from lHaves fc Latham! II1975D. Geneva 7-color photometry 
(U,B1,B , B2,V1,V,G) fromlRufenerl (I1976I). UBVRIJK photome try from lJohnson et all II1966D. mid-IR photome - 
try from Uavawardhana et al.N2001[) and lMetc hey_et alJ l|2003() . IRAS photometry from lBeichman et alJ (^ 
and 60 ^m ISO photometry from Habing et 



itclicv^iaj 
jdj(|2001). 



HR 4796A : The S EP includes optical photometry from lHaves fc Lathaml (119751). UB VRIJHK photometry from 



J_urjLel_aD 1199 3J), KL MN and mid-IR photomet ry from jF^^r^o-Acgste et alJ 1 1998|) , N-band photometry from 
MetehevejT^ilTT2003l) . IRAS photometry from iBeichmanetalT l|1988l) . and submillimeter photometry from 



Jura et alJ <|199 



IRAS 14050- 4109: The SEP includ es optical photometry from IGregorio-Hetem fc Heteml l|2002D . N-band pho- 
tometry from lMetchev et alJ l|2003|) . and IRAS photometry from lBeiclmiaiiet^lTl(T9'88|) . 

HP 16 3296: The SEP includes op tical p hotometry fromlJalnar et alJ l(T995l) . iHauot fc Schrolll II1974D. iThe et~ai] 
( 1985f), and iMalfait et all l|1998l) . a nd lOudmaiier et all ( 120011). UBVRI ph otome try from IHillenbrand et alJ 



1 1.9921) . JHK LMN photometry fro m iBerrilli et alJ (11991. lEiroa et al.l l|2001|) . and lAllenl (I19731). mid-IR pho - 
tometry from lThi et al.l l|200l|) and Uavawardhana et alJ (I2001D. IR AS photo metry from lBeichman et alJ (|l988), 
and (sub)millimeter photometry from iMannings fc Emerson! <|1994) and Oi (200J). 

HP 1 79218: The SEP include s optical photometry from lHaves fc Lathaml (I1975D. UPXYZVS photometry 
fromlEimontas fc Sudziusl lll998tl. TJ BVRI JHKLMN photome try fromlMiroshnichenko et alJ l|1999ft . IMalfait et all 
II1998D. andlLawrence et all (|199(1) . IRAS photometry fro mlBeichman et all II1988D. submillimeter photometry 
from lMeeus et all l|2001|) . and millimeter photometry from Mannings fc Sargentl (|2000fl . 

WW Vul: The SEP includes UPXYZVS photometry fr omlEimontas fc Sudziusl llT99l. IIBVR photom etry from 
iHerbst fc Shevchenkd 1119991) JHK L photometry from lEiroa et all (120011) andlGlass fc Penstonl l|1974D . mid-IR 
phot ometry from lThi et al] ll200i|) . IRAS photometry from IBeichman et al.l l|1988|) . and millimeter photometry 
from Natta et all (I19971I2001D. 



HP 184761: The SEP includ es optical photometry from lHaves fc Lathaml (I1975D. UBVRIJHK photometry from 
IMiroshnichenko et al.l (|1999f) . and IRAS photometry from lBeichman et al.l l|1988|) . 
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HD 2168 03: The SEP i nclud es U PXYZV S phot ometry from lEimontas fc Sudziul (I1998T>. UBVRI phot ome- 
try from I.Tohnson et al.1 l)1966ft and iBessell <ll990lK mid-IR ph otometry from iFaiardo-Acosta et all <ll999ft and 
iMetchev et 811 ( 2003) . IRAS photometry fromlBeichman et all (11988ft . 60 /im ISO photometry from lHabing et alJ 
1)200 If) , and millimeter photometry from lWen"rfcnml7^rirtTi(j^9'^ . 
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